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STATISTICAL METHOD FOR ASSESSING AUTONOMIC BALANCE 

CROSS REFERNCE TO RELATED APPLICATION(S) 
This application is a continuation of U.S. Patent Application No. 10/436,876, 
filed on May 12, 2003, the specification of which is incorporated herein by reference. 

FIELD OF THE INVENTION 
This invention pertains to cardiac rhythm management devices such as 
pacemakers and implantable monitoring devices. 

BACKGROUND 

Heart rate variability (HRV) refers to the changes in the length of time between 
consecutive heart beats during sinus rhythm and is primarily due to the interaction 
between the sympathetic and parasympathetic arms of the autonomic nervous system. 
Measurement and analysis of heart rate variability is thus a useful and non-invasive tool 
for assessing the status of the autonomic nervous system. 

A heart beat is usually measured as the time from the peak of one R wave to the 

peak of the next, referred to as an RR interval. The variability of normal RR intervals 

(i.e., during sinus rhythm) can be determined and analyzed in several different ways in 

either the time domain or the frequency domain. Time domain measurements involve 

the computation of a statistic based upon the individual RR intervals making up an RR 

time series such as the standard deviation of the RR intervals in the series. Frequency 

domain analysis, on the other hand, employs methods such as the Fast Fourier 

Transform (FFT) or autoregressive analysis to analyze the frequency spectrum of the 

variability in the RR intervals. This latter type of analysis has proven to be particularly 

valuable in assessing the relative activities of the sympathetic and parasympathetic 

nervous systems in a subject. Such assessment of the state of autonomic balance would 

be a useful function for implantable cardiac rhythm management devices such as 

pacemakers and implantable cardioverter/defibrillators to perform as it could be used to 

modify the manner in which therapy is delivered by the device or to predict the 
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occurrence of arrhythmias. Frequency domain analysis of heart rate variability, 
however, requires computational and data storage capabilities that may not be practical 
in present-day implantable devices. 

SUMMARY 

The present invention is a method for assessing the autonomic balance of a 
subject by estimating the ratio of certain frequency components in an RR time series 
based upon statistics computed from the RR intervals making up the series. The method 
is especially suitable for use by an implantable device since the required statistics may 
be computed without the processing overhead and data storage capability associated 
with frequency domain analysis. The statistics used to estimate the frequency 
component ratio may be computed by cumulatively summing particular functions of RR 
interval measurements or by a histogram technique in which the relative frequencies of 
RR interval values or functions thereof are calculated from running counts of the RR 
interval measurements whose values are within specified ranges. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a block diagram of an exemplary cardiac rhythm management device 

for practicing the present invention. 

Figs. 2A through 2D illustrate an exemplary spectrum of an RR time series and 

the frequency responses represented by statistical surrogates. 

DETAILED DESCRIPTION 
As noted above, heart rate variability refers to the variability of the time intervals 
between successive heart beats during a sinus rhythm. Spectral analysis of heart rate 
variability involves decomposing a signal representing successive beat-to-beat intervals 
into separate components representing the amplitude of the signal at different oscillation 
frequencies. It has been found that the amount of signal power in a low frequency (LF) 
band ranging from 0.04 to 0.15 Hz is influenced by the levels of activity of both the 
sympathetic and parasympathetic nervous systems, while the amount of signal power in 
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a high frequency band (HF) ranging from 0.15 to 0.40 Hz is primarily a function of 
parasympathetic activity. The ratio of the signal powers, designated as the LF/HF ratio, 
is thus a good indicator of the state of autonomic balance, with a high LF/HF ratio 
indicating increased sympathetic activity. If an implantable medical device could 
5 monitor the LF/HF ratio, the device could log a clinically significant event when the 
ratio exceeds a specified threshold value, as well as possibly automatically altering its 
mode of operation (e.g., initiating different therapies or performing more 
computationally intensive data analysis for arrhythmia prediction). 

A series of RR interval values can be regarded as a discrete signal indexed by 

10 heartbeat such that each value of the signal represents an RR interval for a particular 
heartbeat. In order to properly analyze the frequency content of heart rate variability, 
however, the RR time series should be resampled at a specified sampling frequency in 
order to equalize the time intervals between interval values and thus convert the time 
series into a discrete time signal, where the sampling frequency is selected to meet the 

15 Nyquist criterion with respect to the frequencies of interest. Spectral analysis of such an 
RR interval signal can then be performed directly in the frequency domain using 
discrete Fourier transform or autoregression techniques in order to compute the LF/HF 
ratio. A time-domain technique for determining the high and low frequency 
components of the signal could also be used in which the RR interval signal is input to 

20 low band and high band digital filters and signals proportional to the power of the RR 
interval signal in each of the low frequency and high frequency bands are derived so that 
the LF/HF ratio may be computed. Both frequency domain and time domain analysis 
performed in this manner are computationally intensive, however, and require the 
storage of large amounts of RR interval data. Such methods may therefore not be 

25 practical in a typical implantable medical device which is a small battery-powered 

device with limited processing power. As described below, statistical techniques that do 
not involve such processing overhead may be used to generate surrogate parameters 
from which the LF/HF ratio may be calculated. 
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1 . Exemplary implantable device description 

Cardiac rhythm management devices are implantable devices that provide 
electrical stimulation to selected chambers of the heart in order to treat disorders of 
cardiac rhythm. A pacemaker, for example, is a cardiac rhythm management device that 
paces the heart with timed pacing pulses. The most common condition for which 
pacemakers are used is in the treatment of bradycardia, where the ventricular rate is too 
slow. Cardiac rhythm management devices may also treat tachyarrhythmias, where the 
heart rate is too fast, by anti-tachycardia pacing and/or delivery of defibrillation shocks. 
Such devices are usually implanted subcutaneously on the patient's chest and connected 
to electrodes by leads threaded through the vessels of the upper venous system into the 
heart. An electrode can be incorporated into a sensing channel that generates an 
electrogram signal representing cardiac electrical activity at the electrode site and/or 
incorporated into a pacing or shocking channel for delivering pacing or shock pulses to 
the site. 

A block diagram of an implantable cardiac rhythm management device is shown 
in Fig. 1. The controller of the device is made up of a microprocessor 10 
communicating with a memory 12 via a bidirectional data bus, where the memory 12 
typically comprises a ROM (read-only memory) for program storage and a RAM 
(random-access memory) for data storage. The controller could be implemented by 
other types of logic circuitry (e.g., discrete components or programmable logic arrays) 
using a state machine type of design, but a microprocessor-based system is preferable. 
As used herein, the programming of a controller should be taken to refer to either 
discrete logic circuitry configured to perform particular functions or to executable code 
stored in memory or other storage medium. The controller is capable of operating the 
device so as to deliver a number of different therapies in response to detected cardiac 
activity. A telemetry interface 80 is also provided for enabling the controller to 
communicate with an external programmer. 

The embodiment shown in Fig. 1 has two sensing/pacing channels, where a 
pacing channel is made up of a pulse generator connected to an electrode while a 
sensing channel is made up of the sense amplifier connected to an electrode. A MOS 

4 



Docket No. 279.564US2 

switch matrix 70 controlled by the microprocessor is used to switch the electrodes from 
the input of a sense amplifier to the output of a pulse generator. The switch matrix 70 
also allows the sensing and pacing channels to be configured by the controller with 
different combinations of the available electrodes. The channels may be configured as 

5 either atrial or ventricular channels. In an example configuration, an atrial 

sensing/pacing channel includes ring electrode 43a and tip electrode 43b of bipolar lead 
43c, sense amplifier 41, pulse generator 42, and a channel interface 40. A ventricular 
sensing/pacing channel includes ring electrode 33a and tip electrode 33b of bipolar lead 
33c, sense amplifier 31, pulse generator 32, and a channel interface 30. The channel 

10 interfaces communicate bi-directionally with a port of microprocessor 10 and may 

include analog-to-digital converters for digitizing sensing signal inputs from the sensing 
amplifiers, registers that can be written to for adjusting the gain and threshold values of 
the sensing amplifiers, and registers for controlling the output of pacing pulses and/or 
changing the pacing pulse amplitude. A shock pulse generator 20 is also interfaced to , 

15 the controller for delivering defibrillation shocks through electrodes selected by the 
switch matrix. In the illustrated embodiment, the device is equipped with bipolar leads 
that include two electrodes which are used for outputting a pacing pulse and/or sensing 
intrinsic activity. Other embodiments may employ unipolar leads with single electrodes 
for sensing and pacing. The switch matrix 70 may configure a channel for unipolar 

20 sensing or pacing by referencing an electrode of a unipolar or bipolar lead with the 
device housing or can 60. 

The controller 10 controls the overall operation of the device in accordance with 
programmed instructions stored in memory. The controller 10 interprets electrogram 
signals from the sensing channels in order to control the delivery of paces in accordance 

25 with a pacing mode and/or deliver shock therapy in response to detection of a 

tachyarrhythmia such as ventricular fibrillation. The sensing circuitry of the device 
generates atrial and ventricular electrogram signals from the voltages sensed by the 
electrodes of a particular channel. An electrogram is analogous to a surface ECG and 
indicates the time course and amplitude of cardiac depolarization that occurs during 

30 either an intrinsic or paced beat. When an electrogram signal in an atrial or sensing 
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channel exceeds a specified threshold, the controller detects an atrial or ventricular 
sense, respectively, which pacing algorithms may employ to trigger or inhibit pacing 
and from which heart rates may be derived by measuring the intervals between senses. 

5 2. Statistical assessment of autonomic balance 

A cardiac rhythm management device such as illustrated in Fig. 1 can be 
programmed to determine heart rate variability by analyzing data received from its 
ventricular sensing channels. The intervals between successive ventricular senses, 
referred to as RR intervals, can be measured for a specified period of time or a specified 

10 number of beats and their variability analyzed. A typical RR time series, for example, 
would be made up of RR intervals over 24 hours or other long-term period. In order to 
derive a signal representing heart rate variability during a sinus rhythm, ectopic 
ventricular beats (i.e., premature ventricular contractions or PVCs) can be detected by 
monitoring whether a P wave precedes each R wave, with the RR intervals before and 

15 after the PVC changed to an interpolated or otherwise filtered value. An RR interval 
signal collected in this manner could be spectrally analyzed by the device in order to 
determine the frequency content in the LF and HF bands by either transforming the 
signal into the frequency domain or by decomposing the signal with bandpass filters. 
Both of these approaches, however, require extensive calculations and the storage of 

20 large amounts of data. A less computationally intensive way of spectrally analyzing an 
RR time series is to employ certain statistical parameters as surrogates for the actual 
specific frequency components. 

Fig. 2A illustrates a log frequency plot of the spectrum an example RR time 
series showing a very low frequency component VLF between DC (i.e., no variability in 

25 the RR intervals) and 0.04 Hz, a low frequency component LF between 0.04 Hz and 
0.15 Hz, and a high frequency component HF between 0.15 and 0.40 Hz. Figs. 2B 
through 2D illustrate how the frequency components of an RR time series would be 
reflected by different statistical parameters. The rMSSD statistic is defined as the 
square root of the mean of the squared successive differences of an RR time series: 
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rMSSD = E{(RRj-RRj.i) 2 } 05 

where E is the expectation or mean value operator, and RRj refers to the ith RR interval 
in the series. The square root step in the calculation can be omitted to give the 

5 [rMSSD] 2 parameter. By averaging the square of the successive interval-to-interval 
difference values in the RR time series, the rMSSD or [rMSSD] 2 statistic maximally 
reflects variations in the RR intervals that occur with each successive interval and 
progressively attenuates variations in the RR intervals that occur at lower frequencies. 
The frequency response represented by the rMSSD statistic shown in Fig. 2B is 

10 therefore greatest at the maximum frequency that can be represented in the time series 
(i.e., approximately one-half of the average heart rate, similar to the Nyquist frequency 
in a time series with regular intervals) and then decreases linearly with decreasing 
frequency so that lower frequency variability in the RR time series is not represented. 
Computation of the rMSSD or [rMSSD] 2 statistic thus captures a frequency range of 

15 heart rate variability which is similar to the HF band illustrated in Fig. 2 A. 

The SD 7 statistic is defined as the standard deviation of the mean values of all 
successive 7-second segments in the RR time series, and [SD 7 ] 2 is the square of that 
standard deviation or variance. By computing mean values of the RR intervals over 7- 
second segments, the [SD 7 ] 2 statistic averages out the variations in the RR intervals that 

20 occur over time intervals equal to or smaller than 7 seconds (i.e., variations at 

frequencies higher than approximately 0.15 Hz), with the variations in the RR intervals 
occurring over longer intervals then being reflected by computation of the variance of 
those mean values. Fig. 2C illustrates the frequency response of the [SD 7 ] 2 statistic 
which shows capture of those frequency components of the RR interval signal from 

25 some frequency above DC (since a variance calculation eliminates the DC component of 
a signal) to approximately 1/7 or 0.15 Hz (i.e., a frequency with a period equal to the 
length of the 7-second segment). The SD 2 s statistic is similarly defined as the standard 
deviation of the mean values of all successive 25-second segments in the RR time 
series, where [SD 2 s] 2 is the square of that standard deviation or variance. Fig. 2D 

30 illustrates the frequency response of the [SD 25 ] 2 statistic which shows capture of those 
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frequency components of the RR interval signal from some frequency above DC to 
approximately 1/25 or 0.04 Hz (i.e., a frequency with a period equal to the length of the 
25-second segment). If the [SD 25 ] 2 statistic is subtracted from the [SD 7 ] 2 statistic, the 
resulting parameter captures a frequency range of heart rate variability which is similar 
to the LF band illustrated in Fig. 2A. 

The rMSSD, [SD 25 ] 2 , and [SD 7 ] 2 statistics may thus serve as surrogates for the 
frequency components of an RR time series. An estimate of the LF/HF ratio may then 
be computed as: 

Estimated LF/HF = K {[SD 7 ] 2 - [SD 25 ] 2 }/ [rMSSD] 2 

where K is a constant. A linear regression analysis may be performed in which the 
estimated values are correlated with the actual spectrum of an RR time series to derive 
the value of K. Alternatively, the estimated LF/HF ratio may be compared with 
appropriately scaled threshold values in order to assess the autonomic balance of a 
subject, which eliminates the need for K in the calculation. As described in more detail 
below, estimating the LF/HF ratio in this manner is much less computationally intensive 
than direct spectral analysis since the statistical surrogates for frequency content can be 
computed by maintaining cumulative sums or counts of functions of the measured RR 
intervals. 

A description of an embodiment of the method for estimating an LF/HF ratio 
and assessing a subject's autonomic balance using statistical surrogates which may be 
implemented by an implantable device is as follows. The device measures RR intervals 
between each pair of successive ventricular senses over a predetermined long-term 
period (e.g., 24 hours) to result in RR interval measurements of an RR time series RRi 
through RR N where N is the total number of RR interval measurements during the 
predetermined long-term period. Mean values M x) through M xL of successive x-second 
segments of the RR time series are computed, where x is a predetermined number (e.g., 
7) and L is the total number of such x-second segments in the RR time series. Mean 
values M y i through M yK . of successive y-second segments of the RR time series also 
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computed, where y is a predetermined number greater than x (e.g., 25) and K is the total 
number of such y-second segments in the RR time series. Variances of the mean values 
Mxt through M x l and of the mean values M y i through M yK are next computed, referred 
to as [SD X ] 2 and [SD y ] 2 ; respectively. A mean value of squared successive differences 
5 between the RR intervals in the RR time series is next computed to give the [rMSSD] 
parameter. The ratio of the low frequency content in the RR time series between 
approximately 1/x Hz and 1/y Hz and higher frequency content in the RR time series, 
referred to as LF/HF, can then be calculated as: 



10 LF/HF - K ([SD X ] 2 - [SD y ] 2 ) / [rMSSD] 2 



where K is a defined constant and LF/HF is taken to be reflective of the subject's 
autonomic balance with appropriate selection of the x and y values. 

In order to compute the statistics, a histogram technique may be employed where 

15 cumulative counts of particular functions of the RR interval measurements are 

maintained. For example, to compute the [SD X ] 2 statistic, an integral number T of 
interval bins Ai through At representing interval values Ii through I T are defined. Each 
computed mean value M x j through M xL of the successive x-second segments of the RR 
time series is then assigned to a corresponding one of the interval bins Ai through A T . 

20 That is, each computed mean value is assigned to the bin representing an interval value 
to which it is closest to thereby maintain a cumulative count of the values taken on by 
the computed mean values. The number of computed mean values assigned to each 
interval bin Ai through A T is counted and divided by L to derive a relative frequency 
FXj for each interval value Ii. [SD X ] 2 may then be computed as: 

25 

[SD X ] 2 = Eao'FXi-CECIDFXO 2 



where the summations are carried out from i = 1 to T. 

The [SD y ] 2 statistic may be similarly computed. An integral number S of 

30 interval bins Bi through B s representing interval values Ii through I s are defined. Each 
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computed mean value M yt through M yK of the successive x-second segments of the RR 
time series is then assigned to a corresponding one of the interval bins Bi through B s . 
The number of computed mean values assigned to each interval bin Bi through B s is 
counted and divided by K to derive a relative frequency FYj for each interval value Ij. 
The [SD y ] 2 may then be computed as: 

[SD y ] 2 = E(I0 2 FYME(Ii)FY0 2 

where the summations are carried out from i = 1 to S. 

A histogram technique may also be used to calculate the [rMSSD] 2 statistic. An 
integral number W of RR interval difference bins Ci through C w representing RR 
interval difference values Di through D w are defined. An interval difference between 
each pair of RR intervals in the RR time series is computed as (RR 2 -RRi) through 
(RRn-RRn-i). Each of the computed interval differences between RR intervals in the 
RR time series is then assigned to a corresponding one of the RR interval difference 
bins Ci through C w > and the number of computed interval differences assigned to each 
RR interval difference bin Ci through C w is divided by N-l to derive a relative 
frequency ¥D { for each interval difference Di. The [rMSSD] 2 is then computed as: 

[rMSSD] 2 = E(D0 2 FDi 

where the summation is carried out from i = 1 to W. 

Alternatively, the statistics may be computed by cumulatively summing 
functions of the RR interval differences. For example, the [SD X ] 2 statistic may be 
computed by cumulatively summing each computed mean value M xi and cumulatively 
summing each computed mean value M xi squared so that [SD X ] 2 may be computed as: 

[SD X ] 2 = (l/L)E(M x i) 2 -((l/L)EM xi ) 2 

where the summations are carried out from i = 1 to L by the cumulative summing 
operations. Similarly, the [SD y ] 2 statistic may be computed by cumulatively summing 
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each computed mean value M yi and cumulatively summing each computed mean value 
Myi squared so that [SD y ] 2 may be computed as: 

[SD y ] 2 = (l/L)E(M yi ) 2 -((l/L)EM yi ) 2 

where the summations are carried out from i = 1 to K by the cumulative summing 
operations. The [rMSSD] 2 statistic may be computed by cumulatively summing each 
computed interval difference value (RIWRRi) squared and then computing [rMSSD] 2 
as: 

[rMSSD] 2 =(l/(N-l))E(RRi +r RRi) 2 

where the summation is carried out from i = 1 to N-l by the cumulative summing 
operation. 

As noted earlier, investigators have generally found that the optimum frequency 
ranges for computation of the LF/HF ratio in order to assess autonomic balance is with 
an LF band between 0.04 Hz and 0. 1 5 Hz and with an HF band between 0. 1 5 Hz and 
0.40 Hz. This would correspond to x and y values in the above description of 7 and 25, 
respectively. Different LF and HF frequency ranges and/or different x and y values, 
however, may be found to more optimally assess autonomic balance in a particular 
patient. Also, a standard long-term period over which to define an RR time series is 24 
hours, but a different long-term period may be more appropriate in certain 
circumstances. In another embodiment of the method, a moving average of estimated 
LF/HF ratios over successive long-term time periods may be computed. 

Although the invention has been described in conjunction with the foregoing 
specific embodiments, many alternatives, variations, and modifications will be apparent 
to those of ordinary skill in the art. Other such alternatives, variations, and 
modifications are intended to fall within the scope of the following appended claims. 
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What is claimed is: 

1 . A method, comprising: 

measuring time intervals between successive heart beats, referred to as RR 
5 intervals, over a predetermined long-term period, where N is the total number of RR 
interval measurements during the predetermined long-term period and the measured RR 
intervals RRi through RR N are referred to as an RR time series; 

computing mean values M x \ through M xL of successive x-second segments of 
the RR time series, where x is a predetermined number and L is the total number of such 
10 x-second segments in the RR time series, and computing a variance of the mean values 
M x i through M xL , referred to as [SD X ] 2 ; 

computing mean values M y i through M yK of successive y-second segments of 
the RR time series, where y is a predetermined number greater than x and K is the total 
number of such y-second segments in the RR time series, and computing a variance of 
15 the mean values M y i through M yK , referred to as [SD y ] ; 

computing a mean value of squared successive differences between the RR 
intervals in the RR time series, referred to as [rMSSD] 2 ; and, 

estimating a ratio of the low frequency content in the RR time series between 
approximately 1/x Hz and 1/y Hz and higher frequency content in the RR time series, 
20 referred to as LF/HF, as: 

LF/HF = K ([SD X ] 2 - [SD y ] 2 ) / [rMSSD] 2 
wherein K is a defined constant. 

2. The method of claim 1 further comprising: 

25 defining an integral number T of RR interval bins Ai through A T representing 

RR interval values Ii through It; 

assigning each computed mean value M x i through M xL of the successive x- 
second segments of the RR time series to a corresponding one of the RR interval bins 
A! through At; 
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counting the number of computed mean values assigned to each RR interval bin 
Ai through A T and dividing each such number by L to derive a relative frequency FX* 
for each interval value I,; and, 

computing [SD X ] 2 as 
5 [SD X ] 2 = Eft) 2 FX; -(£(10 FX;) 2 

where the summations are carried out from i = 1 to T. 



3. The method of claim 1 further comprising: 

defining an integral number S of RR interval bins B] through B s representing 
10 RR interval values Ii through Is; 

assigning each computed mean value M y i through M yK of the successive y- 
second segments of the RR time series to a corresponding one of the RR interval bins Bj 
through B s ; 

counting the number of computed mean values assigned to each RR interval bin 
15 Bi through B s and dividing each such number by K to derive a relative frequency FYj 
for each interval value Ijj and, 
computing [SD y ] 2 as 

[SD y ] 2 = £ (Ij) 2 FYj - (E (Ii) FYj) 2 
where the summations are carried out from i = 1 to S. 

20 

4. The method of claim 1 further comprising: 

defining an integral number W of RR interval difference bins Q through Cw 
representing RR interval difference values Di through Dw; 

computing an interval difference between each pair of RR intervals in the RR 
25 time series as (RR 2 -RRi) through (RR n -RRn-i); 

assigning each of the computed interval differences between RR intervals in the 
RR time series to a corresponding one of the RR interval difference bins Ci through Cw; 

counting the number of computed interval differences assigned to each RR 

interval difference bin Ci through Cw and dividing each such number by N-l to derive a 

30 relative frequency FD; for each interval difference Di; and, 

13 



Docket No. 279.564US2 

computing [rMSSD] 2 as 

[rMSSD] 2 = E(D0 2 FDi 
where the summation is carried out from i = 1 to W. 

5 5. The method of claim 1 further comprising: 

cumulatively summing each computed mean value M xi ; 
cumulatively summing each computed mean value M xi squared; 
computing [SD X ] 2 as: 

[SD X ] 2 = (1/L)E(M X1 ) 2 -((1/L)EM X1 ) 2 
10 where the summations are carried out from i = 1 to L by the cumulative summing 
operations. 

6. The method of claim 1 further comprising: 
cumulatively summing each computed mean value My*; 

15 cumulatively summing each computed mean value Myi squared; 

computing [SD y ] 2 as: 

[SD y ] 2 = (l/L)E(M yi ) 2 -((l/L)EM yi ) 2 
where the summations are carried out from i = 1 to K by the cumulative summing 
operations. 

20 

7. The method of claim 1 further comprising: 

cumulatively summing each computed interval difference value (RRj+i-RRO 
squared; 

computing [rMSSD] 2 as: 
25 [rMSSD] 2 = (1/(N-1)) E (RRj+i-RRj) 2 

where the summation is carried out from i = 1 to N-l by the cumulative summing 
operation. 



30 



8. The method of claim 1 wherein the predetermined long-term period is 24 hours. 
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9. The method of claim 1 wherein the predetermined numbers x and y are 7 and 25, 
respectively, such that the estimated LF/HF represents a ratio of the low frequency 
content in the RR time series between approximately 0.04 Hz and 0.15 Hz and the high 
frequency content in the RR time series between approximately 0.15 Hz and 0.4 Hz. 

5 

10. The method of claim 1 further comprising computing a moving average of 
estimated LF/HF ratios over successive long-term time periods. 

11. A system, comprising: 

10 means for measuring time intervals between successive heart beats, referred to as 

RR intervals, over a predetermined long-term period, where N is the total number of RR 
interval measurements during the predetermined long-term period and the measured RR 
intervals RRj through RR N are referred to as an RR time series; 

means for computing mean values M x i through M xL of successive x-second 
15 segments of the RR time series, where x is a predetermined number and L is the total 
number of such x-second segments in the RR time series, and computing a variance of 
the mean values M x i through M xL , referred to as [SD X ] 2 ; 

means for computing mean values M y] through M yK ; of successive y-second 
segments of the RR time series, where y is a predetermined number greater than x and K 
20 is the total number of such y-second segments in the RR time series, and computing a 
variance of the mean values M y i through M y K, referred to as [SD y ] 2 ; 

means for computing a mean value of squared successive differences between 
the RR intervals in the RR time series, referred to as [rMSSD] 2 ; and, 

means for estimating a ratio of the low frequency content in the RR time series 
25 between approximately 1/x Hz and 1/y Hz and higher frequency content in the RR time 
series, referred to as LF/HF, as: 

LF/HF = K ([SD X ] 2 - [SD y ] 2 ) / [rMSSD] 2 
wherein K is a defined constant. 
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12. The system of claim 1 1 further comprising: 

means for defining an integral number T of RR interval bins A, through A T 
representing RR interval values Ii through It; 

means for assigning each computed mean value M x i through M xL of the 
successive x-second segments of the RR time series to a corresponding one of the I 
interval bins Aj through A T ; 

means for counting the number of computed mean values assigned to each ) 
interval bin Ai through A T and dividing each such number by L to derive a relative 
frequency FXi for each interval value l{\ and, 

means for computing [SD X ] 2 as 

[SD X ] 2 = EaO'FXi-^aOFXO 2 
where the summations are carried out from i = 1 to T. 



1 3 . The system of claim 1 1 further comprising: 

means for defining an integral number S of RR interval bins B t through B s 
representing RR interval values Ii through I s ; 

means for assigning each computed mean value M y] through M yK of the 
successive y-second segments of the RR time series to a corresponding one of the RR 
interval bins B t through B s ; 

means for counting the number of computed mean values assigned to each RR 
interval bin Bi through B s and dividing each such number by K to derive a relative 
frequency FYj for each interval value Ii; and, 

means for computing [SD y ] as 

[SD y ] 2 = E (Ij) 2 FYj - (E (Ij) FYj) 2 
where the summations are carried out from i = 1 to S. 
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14. The system of claim 1 1 further comprising: 

means for defining an integral number W of RR interval difference bins Ci 
through C w representing RR interval difference values Di through D w ; 

means for computing an interval difference between each pair of RR intervals in 
5 the RR time series as (RR 2 -RRi ) through (RR N -RR N - 1 ); 

means for assigning each of the computed interval differences between RR 
intervals in the RR time series to a corresponding one of the RR interval difference bins 
Ci through Cw; 

means for counting the number of computed interval differences assigned to 
10 each RR interval difference bin Ci through C w and dividing each such number by N-l 
to derive a relative frequency FD, for each interval difference Dj; and, 
means for computing [rMSSD] 2 as 
[rMSSD] 2 - E(D0 2 FDj 
where the summation is carried out from i = 1 to W. 

15 

1 5 . The system of claim 1 1 further comprising: 

means for cumulatively summing each computed mean value M x j; 

means for cumulatively summing each computed mean value M xi squared; 

means for computing [SD X ] 2 as: 
20 [SD X ] 2 = (1/L) E (M xi ) 2 - ((1/L)E M xi ) 2 

where the summations are carried out from i = 1 to L by the cumulative summing 
operations. 

16. The system of claim 1 1 further comprising: 

25 means for cumulatively summing each computed mean value M y i; 

means for cumulatively summing each computed mean value M yi squared; 
means for computing [SD y ] as: 

[SD y ] 2 = (l/L)E(M yi ) 2 -((l/L)EMyi) 2 
where the summations are carried out from i = 1 to K by the cumulative summing 
30 operations. 
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1 7. The system of claim 1 1 further comprising: 

means for cumulatively summing each computed interval difference value 
(RRi+i-RRi) squared; 

means for computing [rMSSD] as: 
5 [rMSSD] 2 = (1/(N-1)) E (RRi+i-RRi) 2 

where the summation is carried out from i = 1 to N-l by the cumulative summing 
operation. 

18. The system of claim 1 1 wherein the predetermined long-term period is 24 hours. 

10 

19. The system of claim 1 1 wherein the predetermined numbers x and y are 7 and 
25, respectively, such that the estimated LF/HF represents a ratio of the low frequency 
content in the RR time series between approximately 0.04 Hz and 0.15 Hz and the high 
frequency content in the RR time series between approximately 0.15 Hz and 0.4 Hz. 

15 

20. The system of claim 1 1 further comprising means for computing a moving 
average of estimated LF/HF ratios over successive long-term time periods. 
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STATISTICAL METHOD FOR ASSESSING AUTONOMIC BALANCE 

ABSTRACT 

A computationally efficient method for assessing a subject's autonomic balance 
by measurement of heart rate variability is disclosed which is particularly suitable for 
implementation by an implantable medical device. Statistical surrogates are used to 
represent frequency components of an RR time series. A ratio of the low frequency 
component to the high frequency component may then be estimated to assess the 
subject's autonomic balance. 
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